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ABSTRACT 

GRB 990123 was a long, complex gamma-ray burst accompanied by an extremely 
bright optical flash. We find different constraints on the bulk Lorentz of this burst to 
be consistent with the speculation that the optical light is emission from the reverse 
shock component of the external shock. Motivated by this currently favoured idea, 
we compute the prompt reverse shock emission to be expected for bursts in which 
multi-wavelength observations allow the physical parameters to be constrained. We 
find that for reasonable assumptions about the velocity of source expansion, a strong 
optical flash my w 9 was expected from the reverse shocks, which were usually found 
to be mildly relativistic. The best observational prospects for detecting these prompt 
flashes are highlighted, along with the possible reasons for the absence of optical 
prompt detections in ongoing observations. 
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1 INTRODUCTION 

The discovery of a prompt and extremely bright optical flash 
in GRB 990123 (Akerlof et al. 1999), implying an apparent 
peak (isotropic) optical luminosity of 5 x 10 49 erg s _1 , has 
lead to widespread speculation that the observed radiation 
arose from the reverse shock component of the burst. The 
reverse shock propagates into the adiabatically cooled par- 
ticles of the coasting ejecta, decelerating the shell particles 
and shocking the shell material with an amount of inter- 
nal energy comparable to that of the material shocked by 
the forward shock. The typical temperature in the reverse- 
shocked fluid is, however, considerably lower than that of the 
forward-shocked fluid. Consequently, the typical frequency 
of the synchrotron emission from the reverse shock peaks at 
lower energy. It is believed to account for the bright prompt 
optical emission from GRB 990123. The reverse shock emis- 
sion stops once the entire shell has been shocked and the 
reverse shock reaches the inner edge of the fluid. The ejecta 
cool adiabatically after the reverse shock has passed through 
and settled down into a part of the Blandford & McKee 1976 
(BM) solution that determines the late profile of the decel- 
erating shell and the external medium. 

Given the importance of this prediction, we estimate 
the prompt reverse shock emission expected from bursts 
in which multi-frequency data has been able to constrain 
the burst physical parameters. We find a broad range of 
model parameters for which a strong optical flash is expected 
from most of these bursts. Unfortunately, no observations to 
search for an optical flash were performed for these cases. 
The possible reasons for the absence of detectable optical 
flux in other bursts are highlighted, along with the predic- 
tions that may be useful for designing search strategies for 



the rapid follow up of optical flashes. Unlike the continuous 
forward shock, the hydrodynamic evolution of the reverse- 
shocked ejecta is more fragile. As we will demonstrate, the 
temperature of the reverse-shocked fluid is expected to be 
non-relativistic for most of these bursts and thus the evolu- 
tion of their ejecta deviates from the BM solution. The de- 
tection of optical flashes, or firm upper limits, would play an 
important role in discriminating between cold and hot shell 
evolution. Moreover, the strong dependence of the peak time 
of this optical flash on the bulk Lorentz factor F provides 
a way to estimate this elusive parameter. We discuss how 
different constraints on To for GRB 990123 are consistent 
with optical emission from the reverse shock. We assume 
H = 65 kms" 1 Mpc" 1 , fi M = 0.3, and D. A = 0.7. 



2 THE ROLE OF T 

Relativistic source expansion plays a crucial role in virtually 
all current GRB models (Piran 1999; Meszaros 2001). The 
Lorentz factor is not, however, well determined by observa- 
tions. The lack of apparent photon-photon attenuation up 
to ~ 0.1 GeV implies only a lower limit V ~ 30 (Meszaros, 
Laguna & Rees 1993) , while the observed pulse width evolu- 
tion in the gamma-ray phase eliminates scenarios in which 
T >> 10 3 (Lazzati, Ghisellini & Celotti 2001; Ramirez-Ruiz 
& Fenimore 2000). The initial Lorentz factor is set by the 
baryon loading, that is moc 2 , where mo is the mass of the 
expanding ejecta. This energy must be converted to radia- 
tion in an optically-thin region, as the observed bursts are 
non-thermal. The radius of transparency of the ejecta is 
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where _E is the isotropic equivalent energy generated by the 
central site. The highly variable 7-ray light curves can be 
understood in terms of internal shocks produced by velocity 
variations within the relativistic outflow (Rees & Meszaros 
1994). In an unsteady outflow, if F were to vary by a fac- 
tor of ~ 2 on a timescale ST, then internal shocks would 
develop at a distance R± ~ Y 2 c5T > R T . This is followed 
by the development of a blast wave expanding into the ex- 
ternal medium, and a reverse shock moving back into the 
ejecta. The inertia of the swept-up external matter deceler- 
ates the shell ejecta significantly by the time it reaches the 
deceleration radius (Meszaros & Rees 1993), 

=MM 1/3 * 1016 - ioi7cm - (2) 

Given a certain external baryon density no, the initial 
Lorentz factor then strongly determines where both inter- 
nal and external shocks develop. Changes in To will modify 
the dynamics of the shock deceleration and the manifesta- 
tions of the afterglow emission. 



3 REVERSE SHOCK EMISSION AND r 

It has been predicted that the reverse shock produces a 
prompt optical flash brighter than 15th magnitude with rea- 
sonable energy requirements of no more than a few 10 53 
erg emitted isotropically (Meszaros & Rees 1999; Sari & Pi- 
ran 1999a; hereinafter MR99 and SP99). The forward shock 
emission is continuous, but the reverse shock terminates 
once the shock has crossed the shell and the cooling fre- 
quency has dropped below the observed range. The reverse 
shock contains, at the time it crosses the shell, an amount 
of energy comparable to that in the forward one. However, 
its effective temperature is significantly lower (typically by a 
factor of F). Using the shock jump conditions and assuming 
that the electrons and the magnetic field acquire a fraction 
of the equipartition energy e e and eb respectively, one can 
describe the hydrodynamic and magnetic conditions behind 
the shock. 

The reverse shock synchrotron spectrum is deter- 
mined by the ordering of three break frequencies, the self- 
absorption frequency u a , the cooling frequency v c and the 
characteristic synchrotron frequency u m , which are easily 
calculated by comparing them to those of the forward shock 
(MR99; SP99; Panaitescu & Kumar 2000). The equality 
of energy density across the contact discontinuity suggests 
that the magnetic fields in both regions are of comparable 
strength. 

Assuming that the forward and reverse shocks both 
move with a similar Lorentz factor, the reverse shock syn- 
chrotron frequency is given by 

«m = 5.840 x 10 ls 4-i4 / ,-anifo 2 Io,a(l + ^)~ 1/4 Hz, (3) 

while the cooling frequency v c is equal to that of the forward 
shock. Here we adopt the convention Q = 1Q X Q x for ex- 
pressing the physical parameters, using cgs units. The spec- 
tral power F Vm at the characteristic synchrotron frequency 
is 



F Vm = 4.17 J D 2 - 8 2 4 / ,-2^3n^r , 2 (l + ^) 3/8 Jy. (4) 

The distribution of the injected electrons is assumed to 
be a power law of index — p, above a minimum Lorentz 
factor 7^. For an adiabatic blast wave, the correspond- 
ing spectral flux at a given frequency above v m is F„ w 
Fv m { v l v m)~ ~ while below v m is characterised by a 
synchrotron tail with F„ ~ F Vm {i>/v m ) 1 ■ Similar relations 
to those found for a radiative forward shock hold for the 
reverse shock (Kobayashi 2000; hereinafter K00). 

Unlike the synchrotron spectrum, the afterglow light 
curve at a fixed frequency strongly depends on the hydrody- 
namics of the relativistic shell, which determines the tempo- 
ral evolution of the break frequencies v m and z/ c . The forward 
shock is always highly relativistic and thus is successfully de- 
scribed using the relativistic generalisation of the theory of 
supernova remnants. In contrast, the reverse shock can be 
mildly relativistic. In this regime, the shocked shell is unable 
to heat the ejecta to sufficiently high temperatures and its 
evolution deviates from the BM solution (Kobayashi & Sari 
2000; hereinafter KS00). Shells satisfying 

^O.Ol^A-^r-^n-j/^l, (5) 

are likely to have a Newtonian reverse shockQ, otherwise the 
reverse shock is relativistic and it considerably decelerates 
the ejecta. The width of the shell, A, can be inferred directly 
from the observed burst duration by A = cTd ur /(l + z) 
assuming the shell does not undergo significant spreading 
(Piran 1999). 

If £ > 1, then the reverse shock is in the sub-relativistic 
temperature regime for which there are no known ana- 
lytical solutions. In order to constrain the evolution of F 
in this regime it is common to assume V oc R~ 9 where 
3/2 < g < 7/2 (MR99; KS00). For an adiabatic expan- 
sion, T oc T- 9 /( 1+29 > and so u m oc r -3(8+5 S )/7(l+2 fl ) and 

F Vm OC T-(12+ll 9 )/7(l+2 9 ) spectral flux at a given fre _ 

quency expected from the reverse shock gas drops then 

as T -^i9)lnW9) below Um and r -(7+24 P+ 15 P9 )/14(l+2 g ) 

above. For a typical spectral index p = 2.5, the flux decay 
index varies in a relatively narrow range (« 0.4) between 
limiting values of g. 



4 CONSTRAINTS ON THE r OF GRB 990123 

Despite ongoing observational attempts, the optical flash as- 
sociated with GRB 990123 remains the only event of its kind 
detected to date. Observations of this optical flash appear to 
be in good agreement with early predictions for the reverse 
shock emission (Sari & Piran 1999b). As a result, numer- 
ous studies have been done on this event in which reverse 
shock theory has been applied to burst observations in or- 
der to constrain physical parameters and burst properties, 
including To. Current estimates on the bulk Lorentz factor 
for GRB 990123 stretch over nearly an order of magnitude, 
with values ranging from w 200 (SP99) to « 1200 (Wang, 
Dai & Lu 2000). It should be noted, however, that these 

* It should be remarked that equations (3) and (4) refer to this 
reverse shock regime. Equations for the relativistic case are real- 
tively similar, the biggest discrepancy being that the peak flux is 
inversely proportional to T (see equations (7)-(9) of K00). 
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Figure 1. BATSE and ROTSE light curves for GRB 990123 as 
a function of time from the BATSE trigger. Dashed lines rep- 
resent theoretical predictions for the rise oc 

T 3 P -3/2 and decay 

r p-('2\+72,p)/§Q Q £ an a( ji a b a tic reverse shock light curve, as- 
suming the shell is thin and cooling slowly. We predict that the 
optical flash peaked ss 41 ± 6 seconds after the trigger. 



Figure 2. Collective constraints on T for GRB 990123. These in- 
clude estimates from the burst kinematics (narrow distribution), 
synchrotron spectral decay (wider distribution), prompt emission 
pulse width (filled arrow), and jet modelling (unfilled arrow). We 
find a best fit value of T = 900 ± 100. 



1.5+11 (initial 
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estimates were made before accurate burst parameters for 
GRB 990123 were known, and consequently they include ap- 
proximations and parameters from other GRB afterglows. 

By fitting multi-frequency afterglow light curves, phys- 
ical parameters for 8 GRBs have recently been reported 
(Panaitescu & Kumar 2001b, hereinafter PK01). Best fit 
values presented for GRB 990123 are -Ej.so 
jet energy), 6 = 2.1+^, n ,-3 = l-9±i.l > e =.-2 
€b,-4, = 7.4l 23 9 , and p — 2.281q'q3 with a rough estimate 
for the bulk Lorentz factor of To = 1400 ±700 (Panaitescu & 
Kumar 2001a). Using these physical parameters, we present 
a comprehensive examination of the constraints on To and 
report a best-fit value based on an analysis of these con- 
straints. 

Observational estimates for the time of peak flux en- 
abled a measurement of the initial Lorentz factor with rea- 
sonable accuracy using the physical parameters specific to 
GRB 990123. Assuming the optical flash was the result of 
the reverse shock, the initial bulk Lorentz factor 

(1 + zf 8 



237 Elf n~l /a 



rp-3/S 
i d,l 



(6) 



where Tk is the time of peak flux in the observer frame. Light 
curves for the optical flash and 7-ray emission are shown in 
Figure 1. Dashed lines represent theoretical predictions for 
the rise oc T 3p ~ 3/2 and decay oc r - (21+73f,)/96 of an adia- 
batic reverse shock light curve, assuming the shell is thin, i.e. 
A < (E/(2n m p c 2 rf ) )) 1/s , and cooling slowly (K00). Using 
the recently reported physical parameters, one finds GRB 
990123 to have a marginal thickness, as predicted by K00. 
The observed rise time is, however, in good agreement with 
that of a thin shell ^ T 5 ' 5 (in contrast with « T 1/2 for a 
thick shell). The shape of the light curve is determined by 
the time evolution of the three spectral break frequencies, 
which in turn depend on the hydrodynamical evolution of 
the fireball. In the case of GRB 990123, the typical syn- 
chrotron frequency v m = 1.5 x 10 14 Hz is well below the 
cooling frequency v c = 1.0 x 10 19 Hz, and therefore places 
the burst in a regime with a flux decay governed by the rela- 
tion F v oc T - (21+73p)/9(i . This implies a decay of « T~ 2 for 
the optical flash of GRB 990123. Applying these light curve 



predictions to the prompt optical data and taking observa- 
tional uncertainties as well as burst parameter uncertain- 
ties into account, we predict that the optical flash peaked 
Tpcak ~ 41 ± 6 s after the GRB started. Substitution for 
Tpeak in equation (6) gives To = 770 ± 50. 

Observations of the optical peak brightness enable fur- 
ther accurate constraints on the value of To. The syn- 
chrotron spectrum from relativistic electrons comprises four 
power-law segments, separated by three critical frequencies. 
The prompt optical flash in GRB 990123 is observed at a fre- 
quency that falls well below the cooling frequency, but above 
the typical synchrotron frequency: v a < % < ^obs < v c . 
The synchrotron spectrum for this spectral segment is given 
by F obs = F„ m (v obs /v m ) (p - L " 2 (l + 2 ) 1 /2-p/8 wher e u obs is 
taken to be the ROTSE optical frequency. Assuming the 
optical peak flux F obs observed in GRB 990123 is radiation 
arising from the reverse shock, we find To = 1800lgQQ. 

Although observations of a reverse shock induced op- 
tical peak enable fairly accurate calculations of the bulk 
Lorentz factor, it remains possible to obtain information on 
r in situations when an optical flash has not been detected. 
Consider an internal shock which produces an instantaneous 
burst of isotropic 7-ray emission at a time, t, and radius, Ri, 
in the frame of the central engine. The kinematics of collid- 
ing shells implies that although photons are emitted simul- 
taneously, the curvature of the emitting shell spreads the ar- 
rival time of the emission over a period of AT P , thereby pro- 
ducing the observed width in individual pulses. The delay in 
arrival time between on-axis photons and those at 9 ~ 
is a function of the radius of emission and the Lorentz fac- 
tor according to: AT P /(1 + z) = i? ; /(2cr 2 ) where r = T 
in the early phase of the expansion (Fenimore, Madras & 
Nayakshin 1996). In order to allow photons to escape, Ri 
must be larger than the radius of transparency R T . This im- 
poses a lower limit on the initial bulk Lorentz factor such 
that: T > (R T /2cAT p ) 1/2 (l + z) 1/2 . We determine AT P for 
GRB 990123 by measuring the average pulse width through 
autocorrelation methods based on those described in Feni- 
more, Ramirez-Ruiz & Wu (1999) and find AT P « 0.45 s. 
Using the burst parameters to estimate R T from equation 
(1) and applying the inequality relation defined above, we 
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find a lower limit of To > 200. Figure 2 displays the collec- 
tive constraints on To for GRB 990123. An additional lower 
limit of To > 300 constraint from afterglow modelling is also 
included (Panaitescu & Kumar 2002). The combination of 
these constraints leads to an average bulk Lorentz factor for 
GRB 990123 of F = 900 ± 100 and thus R T « 1.3 x 10 14 . 



5 GRB 991023: A UNIQUE OPTICAL FLASH? 

Throughout the two years since its discovery, the optical 
flash associated with GRB 990123 has provoked continu- 
ous study in order to determine whether it was consistent 
with current reverse shock theory. More importantly, if it 
can be proven that this optical flash was the result of a typ- 
ical reverse shock event, then we need to understand why 
it is the only event detected to date. Important informa- 
tion may be gained by estimating the prompt reverse shock 
emission expected from bursts with known physical param- 
eters. In our analysis, we use all 7 bursts from PK01 with 
secure redshifts: 990123; 990510; 991208; 991216; 000301c; 
000926; 010222. Reverse shock analytic light curves were 
derived for each of the bursts to produce a set of predicted 
optical flash spectra given various values of the bulk Lorentz 
factor. Figure 3 displays the reverse shock peak spectra pre- 
dicted for the sub-relativistic regime (Fo = 10 2 , upper pan- 
nel) and the relativistic regime (F = 10 3 , lower panel). The 
mildy-relativistic reverse shock spectrum for GRB 990123 
is included in both panels assuming our best fit value of 
To — 900. The visible ROTSE spectral range is displayed in 
the upper panel, while the multi-frequency bands of Swift 
are shown in the lower panel. It is critical to note that 
the only optical flash which has been successfully detected 
(GRB990123, solid thick line) has one of the lowest predicted 
optical fluxes. If it is assumed that the six GRBs did in fact 
produce optical flashes resulting from reverse shocks, then it 
is clear that their non-detection was not due to comparably 
low peak fluxes. Unfortunately, no early observations were 
performed for this set of bursts. Moreover, the predicted re- 
verse shock emission is too faint in X-rays to be detected 
along with the prompt 7-ray emission. 

Akerlof et al. (2000) reported non-detections of optical 
flashes for six long duration GRBs with localisation errors 
of 1 deg 2 or better (see also Kehoe et al. 2001). One pos- 
sible explanation for these non-detections, apart from both 
positional uncertainties and lack of deep imaging^ could 
be related to the time at which the optical flash reached 
its maximum brightness. It is clear that the time at which 
the reverse shock emission peaks plays an important role 
in the ability to detect a prompt optical flash. Current in- 
struments aimed at detecting flashes are limited by their 
response times, often « 10 seconds after the initial trigger 
(Akerlof et al. 1999; Paczynski 2001). Consequently, there 
exists the possibility that a significant fraction of optical 
flashes peak before images were taken. The peak time can 



t K00 found that the ROTSE limits in at least two of these bursts 
(GRB 981121 and GRB 981223) docs not give strong constraints 
and the lack of detections can easily be explained either by in- 
voking a lower no or if To slighlty deviates from that of GRB 
990123. 
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Figure 3. Estimates of reverse shock peak emission for bursts 
with known physical parameters. The two values of To por- 
tray the two relativistic regimes. The reverse shock spectrum for 
GRB 990123 assumes our best fit value of To = 900. Observed 
SAX/WFC flux for the 990123 7-ray trigger appears in the upper 
panel. ROTSE spectral window and Swift multi-frequency bands 
are shown in the upper and lower panels respectively. 



be estimated with T pea k = max[Td, A/c]. The time of the 
peak in the the reverse shock emission is shown as a func- 
tion of Fo in Fig. 4, along with the lower limits on the bulk 
Lorentz factor reported by Panaitescu & Kumar (2002). We 
find that for reasonable values values of To ~ 10 2 — 10 3 , 50 
per cent of the optical flashes would have peaked before GRB 
990123 and the rest within 3 minutes of their initial trigger. 
Another possible explanation for these non-detections in- 
volves the hydrodynamical evolution of the reverse shock. If 
the reverse shock is only mildly relativistic (£ « 1), then it 
cannot heat the ejecta to sufficiently high temperatures and 
so the Blandford-McKee solution is not valid (KS00). Since 
the flux brightness and photon energy are both directly pro- 
portional to r, determining whether the burst is relativis- 
tic constrains the flux decay. Therefore, to understand and 
detect GRB flashes, we should examine whether events go 
undetected due to the behaviour of the Lorentz factor decay. 
Using PK01 physical parameters and equation (5), we cal- 
culate £ for GRB 990123 and six additional bursts for which 
there was no detected optical flash. Figure 4 (upper panel) 
displays the dependence of £ on To for each burst. Unsur- 
prisingly, we find that the reverse shock in GRB 990123 is 
mildly relativistic with £ = 1.0 ± 0.1, which supports the 
early prediction that £ = 0.7 made by KS00, who stated 
that the value could in fact be larger. The reverse shock in 
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Figure 4. Estimates of £ and T pea ^ for GRBs with available 
multi-frequency data. For reasonable values of To ~ 10 2 — 10 3 , 
we find that: (i) the reverse shock in most cases was expected to 
be mildly relativistic with £ fs 1 (upper panel); (ii) 50 per cent of 
these bursts reached peak emission before GRB 990123 with the 
rest within 3 minutes of their initial trigger (lower panel). 



6 OBSERVATIONAL PROSPECTS AND 
CONCLUSIONS 

An obvious question becomes: how is it that the optical flash 
associated with GRB 990123 remains the only event of its 
kind detected to date? Observations of both faint and heav- 
ily dust-enshrouded optical flashes could be missing due to 
the capability of current instruments. With the launch of 
Swift, the facilities necessary to detect reverse shock emis- 
sion will be made readily available. Within 20-70 seconds 
of the initial trigger, the UV and optical telescope (UVOT) 
will begin collecting images of the burst down to B=24. We 
have predicted that optical flashes brighter than my w 9 for 
at least a subset of bursts, with typical peak times of ~ 50 
seconds and peak frequencies of ps 10 15 Hz. However, there 
exist necessary limitations to our approach: we assumed that 
the equipartition parameters are the same across the con- 
tact discontinuity, and we have used only a subset of bursts 
bearing relatively bright optical afterglows. A more detailed 
analysis of the underlying reasons for the non-detections of 
prompt optical afterglows will require the large and unbiased 
sample of Swift GRB afterglows. 

In summary, we show that reverse shocks are likely to 
produce strong optical flashes in a subset of GRBs, and that 
dust obscuration is a possible reason for the reported non- 
detections. If reverse shock emission turns out to be insignif- 
icant, then the explanation for the bright optical flashes will 
surely be even more remarkable and fascinating. 
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most of these bursts was expected to be mildly relativis- 
tic for Fo < 500. The hydrodynamics of the cold shocked 
ejecta is very different from that of the hot ejecta which 
is described by the BM76 solution. Surprisingly, both cases 
predict rather similar light curves (KS00), with decay laws 
that vary in relatively narrow ranges. Hence, it is unlikely 
that the non-detection of optical flashes is linked to the rel- 
ativistic regime of their reverse shock. 

The absence of detectable optical flux accompanying 
a strong 7-ray emission could also be due to dust obscura- 
tion. There is increasing evidence that GRBs with durations 
longer than 2 seconds are associated with sites of massive 
star formation. This environment will strongly attenuate the 
optical afterglow radiation (Waxman & Draine 2000, Vene- 
mans & Blain 2001, Ramirez-Ruiz, Trentham & Blain 2002). 
The observations of GRB late afterglows are consistent with 
this picture: 50-70 per cent of all bursts have no optical coun- 
terpart down to R=24, which implies a minimum absorption 
of Ah « 2 (Lazzati, Covino & Ghisellini 2002; Reichart & 
Yost 2002). 
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